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I. INTRODUCTION 



The invariant theory of trilinear forms over a three-dimensional complex vector space is 
an old subject with a long history, which, as we shall see, appears even longer if we take 
into account certain indirect but highly relevant contributions 2, 21, 2(1 1^. This question 



has been recently revived in the field of Quantum Information Theory as the problem of 
classifying entanglement patterns of three-qutrit states. 

Indeed, since the advent of quantum computation and quantum cryptography, entangle- 
ment has been promoted to a resource that allows quantum physics to perform tasks that are 
classically impossible. Quantum cryptography P, 111] proved that this gap even exists with 
small systems of two entangled qubits. Furthermore, it is expected that the study of higher 
dimensional systems and of multipartite (e.g. 3-partite) states would lead to more applica- 
tions. A seminal example is the so-called 3-qutrit Aharonov-state, which u is so elegant it 



had to be useful" 



12j : Fitzi, Gisin and Maurer [12J found out that the classically impossible 
Byzantine agreement problem Q| can be solved using 3-partite qutrit states. From a more 
fundamental point a view, the Aharonov state led to non-trivial counterexamples of the con- 
jectures on additivity of the relative entropy of entanglement 2i| and of the output purity 



of quantum channels |37[ . Obviously, these results provide a strong motivation for studying 
3-partite qutrit states. Furthermore, interesting families of higher-dimensional states are 

)erfectly suited to address questions concerning local realism and Bell inequalities (see e.g. 

15^ for a study of three-qutrit correlations). 

It is therefore of interest to find some classification scheme for three-qutrit states. A 



possible direction is to look for classes of equiva 
alent up to local unit ary transformations 111 
SLOCC operations) j] 
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ent states, in the sense that they are equiv- 



32j | or local filtering operations (also called 



321 ]. In the case of three qubits, especially the last 



classification proved to yield a lot of insights (the classification up to local unitaries has too 
much parameters left); the reason for that is that in the closure of each generic orbit induced 
by SLOCC operations, there is a unique state (up to local unitary transformations) with 
maximal entanglement [nl l^. 

In (jy], a numerical method converging to such a maximally entangled state has been 
described. It has been experimentally observed that, when applied to a three qutrit state, 
this method converged to a very special normal form. We shall provide a formal proof of 
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this property, and then study in some detail the geometry of those normal forms. Precise 
statements of the results are summarized in the forthcoming section. 



II. RESULTS 

Let V = C 3 and H = V <8 V g> V regarded representation of the group G = 
SX(3,C) x3 . The elements of 7i will be interpreted either as three-qutrit states 

3 

|V>) = A ^,J,k) (1) 

i ,j,k=l 

or as a trilinear forms 

3 

/ = /(*, y, z) = A mXiyjZk (2) 

i,j,k=l 

that is, we identify the basis state with the monomial XiVjZk- If g = (<? , <7 , <r 3 ^) £ G 
is a triple of matrices, we define x\ = J2 p g-pX p , y] = Y. q 9j}vqA = T, r 9kr x r, and the 
coefficients A!^ k by the condition 

A 'ijk x 'iy'j z k = Yl A ijkXiVjZk, (3) 
the action of G on Ti. being defined by 

9 ■ f = Yl A 'ijk x iVj z k (4) 
It has been shown by Vinberg [33] that a generic state can be reduced to the normal form 



A' x ~ ' I I fK\ 
A ijk = Ud ijk H o e ^' fc O ^ 



— e ijk + — 2" 

(where 5ijk is the Kronecker symbol and the completely antisymetric tensor) by an 
appropriate choice of g G G. 
Our first result is 

Theorem II. 1 When applied to a generic 3-qutrit state (QJ) the numerical algorithm of 
converges to a state which is a Vinberg normal form, generically in the same G-orbit as \ip) . 



As proved in 
More precisely, 



31 



the normal form is unique up to local unitary transformations. 
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Theorem II. 2 A generic state has exactly 648 different normal forms. For special states, 
this number can be reduced to 216, 72, 27 or 1. Moreover, the coefficients u, v, w of the 
normal form can be computed algebraically. 

Theorem II. 3 The coefficients of the normal forms are determined, up to a sign, by an 
algebraic equation of degree 1296, which is explicitly solvable by radicals. 

To form this equation, we need some notions of invariant theory. 

A polynomial P(A) in the coefficients is an invariant of the action of G on Ti if 
P(A') = P{A) for all g G G. These invariants form a graded algebra R (any invariant P is 
a sum of homogeneous invariants) and the first issue is to determine the dimension of the 
space Rd of homogeneous invariant of degree d. The Hilbert series 

h(t) = dim R d t d (6) 

is known j3] 

k[t) = (l-t6)(l_ t 9) (1 _ t 12) ( ? ) 

and in fact, one can prove that R is a polynomial algebra generated by three algebraically 
independent invariants of respective degree 6, 9 and 12. 

The modern way to prove this result is due to Vinberg, who obtained it from his notion 
of Weyl group of a graded Lie algebra, applied to a Z 3 -grading of the exceptional Lie algebra 
E 6 |^|. 

In section HTT1 we shall explain how it can be deduced from the work of Chanler [5j]. We 



prove that certain invariants Iq, Ig and l\2 introduced in Ref. ^ are indeed algebraic gen- 
erators of R and explain how to compute them from the numerical values of the coefficients 
A-ijk, by expressing them in terms of transvectants, that is, by means of certain differential 
polynomials in the form /, rather than in terms of the classical symbolic notation. Given 
the values of the invariants for a particular state, we show how to form and solve the system 
of algebraic equations determining the coefficients, u, v, w of the normal form. 

Let a = Iq, b = Iu and c = I\$ (a certain polynomial in the fundamental invariants). 
Then, the symmetric functions of u 3 , v 3 and w 3 

■0 = u 3 + v 3 + w 3 , x = u 3 v 3 + u 3 w 3 + v 3 w 3 , X = 21Qu 3 v 3 w 3 (8) 
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satisfy 

ip 2 - 12x - a = 

< ^ + Xip-b = (9) 
^6 _ | A ^,3 _ i A 2 _ c = o 

Theorem II. 4 The system (QJ) /ias generically 1296 solutions (u,v,w), which can be ob- 
tained by solving a chain of algebraic equations of degree at most 4. Only 648 of them give 
the correct sign for Ig . The number of solutions ( with the correct sign for Ig ) can be reduced 
only to 216, 72, 27 or 1. Moreover, the isotropy groups of these degenerate orbits can be 
determined, and the configuration of the points (it, v, w) in C 3 can be interpreted in terms 
of the geometry of regular complex polyhedra. 

The details are given in Section IYIII 



III. THE FUNDAMENTAL INVARIANTS 



In this section, we describe the fundamental invariants, as well as the other concomitants 
obtained by Chanler p], in a form suitable for calculations, in particular for their numerical 
evaluation. 

As already mentioned, we shall identify a three qutrit state G TC with a trilinear form 

/(x,y,z)= ^ AijkXiVjZk (10) 

l<i,j,k<3 

in three ternary variables. To construct its fundamental invariants, we shall need the notion 
of a transvectant, which is defined by means of Cayley's Omega process (see, e.g., Ref. 

Let fx, /2 and f% be three forms in a ternary variable x = (x 1? x 2 , x 3 ). Their tensor 
product fi Cg) f 2 ® /3 is identified with the polynomial f 1 {K^)f 2 (x^)f 3 (x^) in the three 
independent ternary variables x 1 - 1 ^, x 1 - 2 - 1 and x^ 3 -*. We use the "trace" notation of Olver 
to denote the multiplication map fi®f2®fz^ f 1/2/3, that is, 

tr/ 1 (x( 1 ))/ 2 (x( 2 ))/ 3 (x( 3 )) = / 1 (x)/ 2 (x)/ 3 (x) (11) 
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Cayley's operator fi x is the differential operator 
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8x^ 


a*< 3 > 
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a 






&<»> 



(12) 



Now, we consider three independent ternary variables x, y and z together with the associated 
dual (contravariant) variables £ = (£1,62, £3), V = (^1,^2,^3), C = ((1^2,(3) (that is, £, is 
the linear form on the x space such that £i(x_j) = Sij). 

A concomitant of / is, by definition, a polynomial F in the A^ x, y, z, £, 77, £, such 
that if # = (<?i, g 2 , 53) G SX(3,C) 3 , then, with A', x' etc. as above, 



F(A; x, y , z ; £', 77', C') = x, y, z; £, 77, £). 



(13) 



The algebra of concomitants admits only one generator of degree 1 in the Aijk, which 
is the form / itself. Other concomitants can be deduced from / and the three absolute 
invariants P a = ^£i£;, Pp = J^VjVj an d Pj = J2(k z k, using transvectants. If Fx, F 2 and 
F 3 are three 6-tuple forms in the independent ternary variables x, y, z, £, rj and (, one 
defines for any (m, n 2 , ns) X (mi, 777.2, 7713) £ N 3 x N 3 the multiple transvectant of Fi, F 2 and 
^3 by 



(F 1 ,F 2 ,F s ) nin2n3 =tr ^m^n^n™ 1 ^ 2 ^ 3 



nF(x W ,y W ,z«;£«,77W,C ( ^ 



(14) 



For convenience, we will set (F 1} F 2 , F 3 ) nin2ns = (F x , F 2 , F 3 )^ 2n3 . The concomitants of 
degree 2 given by Chanler P can be obtained using these operations: 



Qa — {f, f, PpP-) 



\011 



(fJ,P a P ( ) 101 



Q-y — if j /) -fo-f/ 



aio 



002 



The invariant Jg is then 

^6 = gg(Q<*' Qm Qa)ln - gg(Qy3, <5/9,<5^)l01 = gg(<^ 

There is an alternative expression using only the ground form /, 

1 



1152 



(15) 
(16) 
(17) 

(18) 

(19) 
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Now, in degree 3 the covariants B a , Bp and P 7 of Ref. are 

B a = (fJJ) 0U (20) 

Bp =(fJJ) im (21) 

B, = (fJJ) lw . (22) 

The other concomitants found by Chanler can be written in a similar way: 

C aP = \(fJJP,) 1W (23) 

Cp a = \(fJJP a ) no (24) 

C ai = \ifJJP,) 101 (25) 

C ia = \(fJJP a ) W1 (26) 

C„ = \ifJJPT 1 (27) 

C lP = lifJJP,) 011 (28) 

D a = -2(fPp,fP„f) ul (29) 

Dp = 2(/P Q ,/P 7 ,/) m (30) 

D, = -2(fP a JPp,f) in (31) 

E a = (Q a J,P a ) 100 (32) 

E = (Qp,f,P p ) 0W (33) 

P 7 = (Q 7 ,/,P 7 ) 001 (34) 

G a = --(.f^, /P 7 , /) 011 + ^(/P„P 7 , /, /)° n (35) 

C/3 = -|(/P«, fP 7 , /) 101 + ^(/P«P 7 , /, /) 101 (36) 

G 7 = -|(/P„, /P* /) 110 + ^(/P.P/3, /, /) 110 (37) 

P = ±(fP a JPp,fP 7 ) ni (38) 



Here, we have combined the concomitants of degree 0, 1 and 2 into independent concomitants 
of degree 3. Next, we have chosen the scalar factors so that the syzygies given by Chanler 



15] 



hold in the form 



H + E a — Ej + DpPp = 0, (39) 

H+E(3-E a + Dry Pry = 0, (40) 

H + E J -E (3 + D a P a = 0, (41) 

3C a p - B 7 P P = 0, (42) 

SCfla - BryP a = 0, (43) 

3C Q7 - £gP 7 = 0, (44) 

3C 7Q - B p P a = 0, (45) 

3C M - B a Pry = 0, (46) 

3C 7/3 - = 0, (47) 

6G Q - 3Q a f + B a P p P n = 0, (48) 

6^ - 3Qfsf + BpP a Pry = 0, (49) 

6G 7 - 3Q 7 / + B^Pp = 0. (50) 

One can remark that a basis of the space of the concomitants of degree 3 found by Chanler 
can be constructed using only transvections and products from smaller degrees, 

/ ) Qafl Qpf, Q7/; B a , Bf3, Bry, D a , Dp, Dry, E Q . (51) 

The knowledge of these concomitants allows one to construct the invariants Ig and 7x2 

I9 = —(E a ,Ep,Ei3)\ll (52) 

In = ^^{B a f,B a f,B a fY n . (53) 

These expressions, which can be easily implemented in any computer algebra system, will 
prove convenient to compute the specializations discussed in the sequel. 

IV. NORMAL FORM AND INVARIANTS 

It will now be shown that a generic state can be reduced to the normal form 

. w — V w + v 

Aijk — udijk H — ejjfc H — \6ijk\ (54) 
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where is the alternating tensor, or, otherwise said, that the generic trilinear form 
/(x, y, z) is equivalent to some 

N uvw (x, y, z) = uixiyiZi + x 2 y 2 z 2 + x 3 y 3 z 3 ) 

+v(x 1 y 3 z 2 + x 2 y x z 3 + x 3 y 2 z x ) (55) 
+w(xiy 2 z 3 + x 2 y 3 z x + x 3 yxz 2 ) 

For such a state, the local density operators are all proportional to the identity. This property 
is automatically satistfied by the limiting state obtained from the numerical method of Ref. 
3l| , and implies maximal entanglement as well. Since this algorithm amounts to an infinite 
sequence of invertible local filtering operations, the genericity of Vinberg's normal form, 
together with the previously mentioned properties, implies convergence to a Vinberg normal 
form for a generic input state, that is, our theorem III. 11 

This normal form is in general not unique, and the relations between the various N uvw in 
a given orbit is an interesting question, which will be addressed in the sequel. 

Although the validity of this normal form follows from Vinberg's theory js^j], it can also 
be proved in other ways, some of them being particularly instructive. We shall detail one 
of these possibilities, which will give us the opportunity to introduce some important poly- 
nomials, playing a role in the algebraic calculation of the normal form and in the geometric 
discussion of the orbits. 

The shortest possibility, although not the most elementary, relies on the results of Ref. 
0| , and starts with computing the invariants of N uvw . We then use a few results of algebraic 
geometry, which can be found in 3^|. Let us denote by Ck = Ck{u,v,w) (k = 6,9, 12) the 
values of the on N uvw . Direct calculation gives, denoting by m pqr the monomial symmetric 
functions of u,v,w (sum of all distinct permutations of the monomial u p v q w r ) 

C & = 771(6) - 10771(3,3) , (56) 

C 9 = (u 3 - v 3 ) (u 3 - w 3 ) (v 3 - w 3 ) , (57) 
Cxi = ro(i2) + 4to( 9)3 ) + 6m (6j6) + 228m {6i 3 i3) . (58) 

It is easily checked by direct calculation that the Jacobian of these three functions is nonzero 
for generic values of (u, v, w). Actually, its zero set consists of twelve planes, whose geometric 
significance will be discussed below. 
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Let us denote by ip : TC ^ Ie ' l9 ' Il2 \ (£ 3 ? fog ma p sending a trilinear form to its three 
invariants, so that (Cq,Cq,Cu) = ip(N uvw ). Let S = {N uvw \(u, v , w) 6 C 3 } be the three 
dimensional space of normal forms. The nonvanishing of the Jacobian proves that (p induces 
a dominant mapping from S to C 3 (that is, the direct image of any non-empty open subset 
of S contains a non-empty open subset of C 3 ). Note that the independence of C 6 ,C 9 ,Ci 2 
implies the independence of Iq, Ig, I\2- Now, Chanler p has shown that 16,-^9, 112 separate 
the orbits in general position. This proves that the field of rational invariants of G is fr eely 
generated by Ig, I 9 ,ii 2 (|35j. Lemma 2.1). As a consequence, (p is a rational quotient ([35], 
section 2.4) for the action of G on 7i (actually, this also implies that ip is a categorical 
quotient, by 3^|, Proposition 2.5 and Theorem 4.12, using that tp\$ is surjective, whence 
also tp). 

There exists a non-empty open subset Yq of C 3 such that the fiber of ip over each of its 
points is the closure of an orbit ( 35], Proposition 2.5). Let then Uq = y? _1 (Y"o)- This set 
cuts S since <p\s is dominant. Let U\ be the union of all orbits having maximal dimension 
(a nonempty open set, the function dimension of the orbit being lower semi-continuous). It 
is easy to see that U\ intersects S (for instance at u = 1, v = 1, w — —1, whose orbit has 
dimension 24 = dimG, as may be checked by direct calculation). Let Sq = U\ fl S, a dense 
open subset of S. The set ip~ 1 {p(So) thus contains a dense open subset U2 of H. One then 
checks that UoClUiD U2 (a dense open subset, as an intersection of dense open subsets of an 
irreducible space) is contained in GS. This proves GS = 7i, that is, the normal form N uvw 
is generic. 

Let us remark that the ahove dlSCUSsi o n also proves, thanhs to Igusa's theoten, (Q, The- 
orem 4.12) that C [H] G = C [Iq, 1$, 1 12], that is, the algebra of invariants is freely generated 
by Chanler's invariants. 

Is is also possible to give a direct proof of the normal form by using the same technique 
as in Ref. |5(. Chanler's method rely on the geometry of plane cubics, which will play a 
prominent role in the sequel. 
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V. THE FUNDAMENTAL CUBICS 

The trilinear form /(x, y, z) can be encoded in three ways by a 3 x 3 matrix of linear 
forms M x (x), M y (y) and M z (z), defined by 

/(x, y, z) = VM.ixjz = 'xM y (y)z = *xM z (z)y (59) 

and the classification of trilinear forms amounts to the classification of one of these matrices, 
say M x (x), up to left and right multiplication by elements of SL(3, C ) and action of SL(3, C ) 
on the variable x. 

The most immediate covariants of / are the determinants of these matrices 

X(x) =deti\4(x) = iS Q , (60) 
Y(y) =detM y (y) = lB p , (61) 
Z(z) = detM 2 (z) = ±B y . (62) 

These are ternary cubic forms, and for generic / the equations X(x) = 0, etc. will define 
non singular cubics (elliptic curves) in P 2 . It is shown in Ref. |29| that whenever one of 
these curves is elliptic, so are the other two ones, and moreover, all three are projectively 
equivalent. Actually, one can check by direct calculation that they have the same invariants. 
When / = N umu , these three cubics have even the same equation and are in the Hesse 
canonical form [6J] 

X(x) = -4>(xl + x\ + x\) + ip x x x 2 xz = Y (x) = Z(x) (63) 



where we introduced, following the notation of Ref. |2ll |. 

<p = uvw , ip = u 3 + v 3 + w 3 . (64) 

The Aronhold invariants of the cubics ()63|) are given by 

Q 4 S = -<f){4j 3 + (60) 3 ) (65) 

6 6 T= (60) 6 + 2O(60 3 )^ 3 - 8ip 6 . (66) 

These are of course invariants of /. We recognize that 6 4 S* = — C12, and we introduce an 
invariant lis such that Ci& = Iis(N uvw ) = 6 6 T. The three cubics have the same discriminant 
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6AS 3 + T 2 , known to be proportional to the hyperdeterminant of / (see Refs. 22 1), which 
we normalize as 

A = 27(64S 3 + T 2 ) . (67) 
Then A = C'\ 2 , where C'vi is the product of twelve linear forms 

C'\2 = uvw(u + v + w)(eu + v + w)(u + ev + w) 
x (e 2 u + ev + w) (u + e 2 v + w) 
x (eu + ev + w) (e 2 u + v + w) 
x(eu + e 2 v + w)(e 2 u + e 2 v + w) (68) 

where e = e 2l7f//3 , so that C\i = is the equation in P 2 of the twelve lines containing 3 by 3 
the nine inflection points of the pencil of cubics 

u 3 + v 3 + w 3 + 6m uvw = (69) 

obtained from X, Y, Z by treating the original variables as parameters. We note also that 
the Jacobian of Ce,Cg, C\2 is proportional to C'\ 2 - 

VI. SYMMETRIES OF THE NORMAL FORMS 



In this section, we will prove Theorem 111.21 That is, a generic / has 648 different normal 
forms ( the cases where this number is reduced will be studied in Section fVIIj) . 

To prove the theorem, we remark that the Hilbert series (J7J) is also the one of the ring of 
invariants of G25, the group number 25 in the classification of irreducible complex reflection 
groups of Shephard m d Todd Q. This group, „h,ch we wffl denote for short by K, has 
order 648. It is one of the groups considered by Maschke |2J| in his determination of the 
invariants of the symmetry group of the 27 lines of a general cubic surface in P 3 (a group 
with 51840 elements, which is related to the exceptional root system E e ). To define K, we 
first have to introduce Maschke's group H, a group of order 1296, which is generated by the 
matrices of the linear transformations on C 3 given in Table 1. 

This group contains in particular the permutation matrices, and simultaneous multipli- 
cation by ±e k , since E 2 = —B. The subgroup K is the one in which odd permutations can 
appear only with a minus sign. It is generated by A, C, D, E. 
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TABLE I: The generators of H 





A 


D 


C 


D 


E 


u' 
v' 
w' 


V 

w 
u 


u 
w 

V 


u 
ev 
e 2 w 


u 

ev 
ew 


-^(u + v + w) 
-^{u + ev + e 2 w) 
-^(u + e 2 v + ew) 



Then, as proved by Maschke, the algebra of invariants of K in C [u, v, w] is precisely 

C [Cq, Cg, C12). 

Hence, we can conclude that K is the symmetry group of the normal forms N uvw . There 
was another, equally natural possibility leading to the same Hilbert series. The symmetry 
group L of the equianharmonic cubic surface E : Zq + z\ + z\ + z\ = acting on the 
homogeneous coordinate ring C [£] has as fundamental invariants the elementary symmetric 
functions of the zf, the first one being by definition, so that the Hilbert series of C [S] L 
coincides with (JJJ). Moreover, L is also of order 648, but it is known that it is not isomorphic 
to K. 

Taking into account the results of Section IIV1 we see that 

S = {N uvw \(u,v jW )eC 3 } (70) 

is what is usually called a Chevalley section of the action of G on TC, with Weyl group K 
(see Ref. (3^, p. 174). This implies that each generic orbits intersects S along a K orbit, 
which in turn implies Theorem 111.21 



VII. THE FORM PROBLEM 

This section contains the proofs of Theorem 111.31 and 111.41 We shall formulate and solve 
what Felix Klein (see Ref. |l6| ) called the "Formenproblem" associated to a finite group 
action. This is the following: given the numerical values of the invariants, compute the 
coordinates of a point of the corresponding orbit. 

In our case, we shall see that the problem of finding the parameters (u, v, w) of the normal 
form of a given generic /, given the values of the invariants, can be reduced to a chain of 
algebraic equations of degree at most 4, hence sovable by radicals. 

Let a = Iq, b = I12 and c = lis (we start with lis, because Gig is a symmetric function of 
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u 3 ,v 3 ,w 3 , and at the end of the calculation, select the solutions which give the correct sign 
for Cg, which is alternating). 

What we have to do is to determine the elementary symmetric functions e\ — ip, = 
X, e 3 = 3 of u 3 , v 3 , w 3 . Let A = 2160 3 . Then, 

^ + X^-b =0 (71) 
^-\\^- l -\ 2 -c =0. (72) 

Eliminating A from these equations, we get a quartic equation for ip 2 

27^j 8 - 186V 4 - 8ap 2 - b 2 = . (73) 

The discriminant (with respect to ip) of this polynomial is proportional to D = b 2 (b 3 — c 2 ) 4 . 
When it is non zero, we get 8 values for e±, each of which determines univocally €2 and 
e 3 . Hence, we obtain eight cubic equations whose roots are the possible values of u 3 ,v 3 ,w 3 . 
This gives 8 sets, whence 8 x 6 = 48 triples, each of which providing generically 27 values 
of (u, v, w), in all 48 x 27 = 1296 triples corresponding to the given values of a, b, c, among 
which exactly 1296/2 = 648 give the correct sign for Ig. The common discriminant of the 
8 cubics is 5 = a 3 — 3ab + 2c. Clearly, when 5 ^ 0, we will have 648 triples. If 5 = 0, one 
can check that the cubics cannot have a triple root, and that no root is zero. Hence, in this 
case, we obtain again 648 triples. 

If D = 0, we can have b 3 = c 2 or b = 0. In the first case, setting b = q 2 ,c = q 3 , the 
equation becomes 

(^ 2 -g) 3 (^ 2 + 2g) 3 = 0. (74) 

In this case, we get only four quartics for if) 2 . If Cg ^ 0, we obtain 216 triples. If Cg = 
and b = a 2 /4, c = — a 3 /8 we obtain again 216 triples which form the centers of the edges 
of a complex polyhedron of type 2{4}3{3}3 in C 3 (see Fig. [TJ, in the notation of Ref. |8j. 
The vertices of this polyhedron are the vertices of two reciprocal Hessian polyhedra (see 
Fig. |2J) and its edges join each vertex of one Hessian polyhedron to the 8 closest vertices 
of the other one. In Fig. the edges of the Hessian polyhedron, which are complex lines, 
are represented by real equilateral triangles, so that the figure can as well be interpreted as 
a 2-dimensional projection of a 6-dimensional Gosset polytope 22i- If Cg = and b = a 2 , 
c = a 3 , we obtain only 72 triples which are the centers of the edges of a Hessian polyhedron 
and the vertices of a complex polytope of type 3{3}3{4}2 (see Fig. EJ). 
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FIG. 2: The Hessian polyhedron 

In the case where b = 0, we have to distinguish between the cases c / and c = 0. If 
c ^ 0, we find 648 triples, whatever the value of a. If c = 0, we obtain 27 triples if a ^ 0, 
and only one if a = 0. 

Indeed, for b = c = 0, the V'-equation reduces to ip 8 = 0, and all the cubics collapse to 
12U 3 — all = 0. For a^Owe obtain precisely 27 triples (u, v, w) which form the vertices of 
a Hessian polyhedron in C 3 (see Ref. p]). 

From the results of Ref. [26j about the arrangement of 12 planes formed by the mirrors 
of the pseudoreflections of K = G25, we can determine the structure of the stabilizers of the 
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FIG. 3: The polyhedron 3{3}3{4}2 

normal forms. The only nontrivial cases are: 

• the orbits with 216 elements, for which the stabilizer is the cyclic group C3; 

• the orbits with 72 elements, for which it is C3 x C3; 

• the Hessian orbits with 27 elements, for which it is the group G4 of the Shephard-Todd 
classification. 

These results can be regarded as a complete description of the moduli space of three 
qutrits states. To see what this means, let us recall some definitions from geometric invariant 
theory. 

It is well known that it in general, the orbits of a group action on an algebraic variety 
cannot be regarded as the points of an algebraic variety. To remedy this situation, one has 
to discard certain degenerate orbits. It is then possible to construct a categorical quotient 
and a moduli space, which describe the geometry of sufficiently generic orbits, respectively 
in the affine and projective situation. 

The categorical quotient Y = TC//G is defined as the affine variety whose affine coordinate 
ring is the ring of polynomial invariants R = C [H] G . The moduli space is the projective 
variety M. = Proj (R) of which R is the homogeneous coordinate ring. It is the quotient of 
the set P(7-^) ss of semi-stable points by the action of G (by definition, a point is semi-stable 
iff at least one of its algebraic invariants is nonzero, see |35{| ) . 
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Now, since in our case the algebra of invariants is a polynomial algebra, we see that the 
categorical quotient is just the affine space C 3 . 

The moduli space is more interesting. The projective variety whose homogeneous coor- 
dinate ring is a polynomial algebra over generators of respective degrees di, . . . , d m is called 
a weighted projective space P(di, . . . , d m ). Hence, by definition, our moduli space M. is the 
weighted projective space P(6, 9, 12) ~ P(2, 3,4). It is known that this space is isomorphic 
to P(l,2,3) |9|, which in turn can be embedded as a sextic surface in P 6 , the so-called del 
Pezzo surface F 6 (see Ref. |l4j)- The del Pezzo surfaces are very interesting objects, known 
to be related to the exceptional root systems (see, e.g., Ref. [20]). 

The above results can then be interpreted as a description of the singularities of Ai, since 
one can view it as the quotient of the projective plane P 2 of the parameters (u : v : w) under 
the projective action of (?25- We have described this quotient as a 648-fold ramified covering 
P 2 — > Ai, and analyzed its ramification locus. 



VIII. CONCLUSION 



A problem of current interest in Quantum Information Theory has been connected to 
various important mathematical works, scattered on a period of more than one century 
from Ref. j^lj in 1889 to Ref. j3| in 2000, in general independent of each other and 
apparently discussing different subjects. Relying on all these works, we have described the 
geometry of the normal forms of semi-stable orbits of three qutrit states under the action of 
SL(3, C ) x3 , the group of local filtering (SLOCC) operations. From a physical point a view, 
our results can be expected to provide a good starting point for studying the richness of 
the entanglement of three qutrits and its differences with that of the simpler qubit systems. 
From a mathematical point a view, we have worked out an interesting example of a problem 
in invariant theory, using both classical algebraic and modern geometric methods, found a 
surprising connection with the geometry of complex polytopes, and applied Klein's vision 
of Galois theory to the explicit solution of an algebraic equation of degree 648. 

Also, this example provides a good illustration of the ideas presented in Refs. (2] and 
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